Unidentified 7-ray sources off the Galactic plane as low- 
mass microquasars? 
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j-jT). Abstract. A subset of the unidentified EGRET 7-ray sources with no active galac- 

i—{ ' tic nucleus or other conspicuous counterpart appears to be concentrated at medium 

latitudes. Their long-term variability and their spatial distribution indicate that they 
are distinct from the more persistent sources associated with the nearby Gould Belt. 
s ! ' They exhibit a large scale height of 1.3 ±0.6 kpc above the Galactic plane. Potential 

counterparts for these sources include microquasars accreting from a low-mass star 
and spewing a continuous jet. Detailed calculations have been performed of the jet 
inverse Compton emission in the radiation fields from the star, the accretion disc, 
and a hot corona. Different jet Lorentz factors, powers, and aspect angles have been 
explored. The up-scattered emission from the corona predominates below 100 MeV 
] whereas the disc and stellar contributions are preponderant at higher energies for 

qq . moderate (~ 15°) and small (~ 1°) aspect angles, respectively. Yet, unlike in the 

high-mass, brighter versions of these systems, the external Compton emission largely 
, fails to produce the luminosities required for 5 to 10 kpc distant EGRET sources. 

Synchrotron-self-Compton emission appears as a promising alternative. 
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1. Variable 7-ray sources in the Galactic disc 

The EGRET telescope has detected 263 7-ray sources above 100 MeV 
(Hartman et al., 1999), half of which have been firmly or plausibly 
associated with flat radio spectrum AGN and a handful of nearby 
radiogalaxies (Sowards-Emmerd et al., 2003 ; Sowards-Emmerd et al., 
2004; Mattox et al., 2001). Six are identified with pulsars. The other 126 
sources remain unidentified because of their poor localization, typically 
within a degree, or their faintness at lower energies. 

Thirty three (33) bright sources along the Galactic plane are associ- 
ated with star-forming regions that harbour many likely 7-ray emitters, 
steady ones such as pulsars and supernova remnants, and more variable 



f Fellow of CONICET 
§ Member of CONICET 



© 2008 Kluwer Academic Publishers. Printed in the Netherlands. 



Grenier_Kaufman_HK04_rev.tex; 2/02/2008; 12:38; p.l 



2 



ones such as pulsar wind nebulae, massive star binaries, and accreting 
X-ray binaries (Romero et al., 1999; Grenier, 2004). Nolan et al. (2003) 
found a group of 17 sources, concentrated in the inner Galaxy, that 
exhibit variability on timescales of weeks to months, hence ruling out 
a supernova remnant or a young pulsar origin. This group is brighter 
than the other low-latitude sources (Bosch- Ramon et al., 2004). Their 
luminosity ranges from 0.8 to 20 x 10 34 {D / 5 kpc) 2 erg s" 1 sr" 1 above 
100 MeV. They may be associated with variable pulsar wind nebulae 
or with lower-energy versions of the colliding winds from a pulsar and 
a massive star, as recently observed at TeV energies by HESS in PSR 
B1259-63 (Kirk et al., 1999). 

Kaufman et al. (2002) have alternatively suggested that they are mi- 
croquasars with persistent jets and high-mass stellar companions. The 
long-term variability may be caused by the precession of the accretion 
disc, henceforth of the jet, that is induced by the gravitational torque of 
the star, or by variations in the accretion rate along an eccentric orbit, 
or by instabilities and inhomogeneities in the stellar wind. Synchro-self- 
Compton (SSC) radiation in the expanding radio blobs ejected during 
sporadic flares shine too briefly to account for EGRET sources (Atoyan 
& Aharonian, 1999). Three microquasars with high-mass companions 
indeed coincide with EGRET sources: LS5039 in 3EG J1824-1514, LSI 
+61 303 in 3EG J0241+6103, and AX J1639.0-4642 in 3EG J1639-4702 
(see Ribo et al, these proceedings). The former and the latter do not 
exhibit variability in the EGRET data (Nolan et al., 2003). The jet 
from LSI +61 303 is known to precess and the observed 7-ray flux is 
variable, though it does not correlate with the radio flares (Kniffen et 
al. (1997); see , however, Massi (2004)). Energetic electrons or pairs 
in the jet can spawn 7-rays by up-scattering photons from the bright 
UV star, the soft X-ray accretion disc, a hard X-ray corona, or their 
own synchrotron radiation. External inverse Compton (EC) emission 
from a cylindrical jet (Romero et al., 2004), as well as SSC emission 
or EC emission from an expanding jet (Bosch- Ramon et al., 2004), can 
both reproduce the observed luminosities. SSC emission predominates 
for large jet magnetic fields, typically > 100 G. EC emission yields a 
large variety of 7-ray spectra, depending on the relative contributions 
of the different radiation fields and on the jet/accretion power ratio. 
The latter should typically exceed 10~ 4 to power an EGRET source. 



2. Variable 7-ray sources off the Galactic plane 

The 93 unidentified sources detected away from the Galactic plane are 
displayed in Figure 1. Their concentration at 3° < \b\ < 30° and in the 
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Figure 1, All-sky plot, in Galactic coordinates, of the unidentified EGRET sources 
at latitudes |6| > 3°. The steadier sources associated with the Gould Belt and the 
variable sources are marked as circles and stars, respectively. 
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Figure 2. Distributions of the S variability indices from Nolan et al. (2003) for the 
firm AGN, Gould Belt, and halo sources at latitudes \b\ > 3°. 



inner half steradian clearly indicates (at a la confidence level) that 70 
to 100 % of them have a Galactic origin, depending on the choice of 
Galactic scale height. Their temporal and spatial characteristics reveal 
a heterogeneous sample. A subset of steadier sources, marked as circles 
in Figure 1, has been associated with the nearby Gould Belt (Gehrels et 
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al., 2000; Grenier, 2000). The others are clearly more variable. Figure 
2 shows their variability index distribution which closely follows that 
of the variable firm AGN sources (Nolan et al., 2003). The two dis- 
tributions are consistent under the Kolmogorov (K) test or the T test 
developed by Eadie et al. (1977), and both are significantly at variance 
with the Belt source distribution (5a and 4a). The average S indices 
of the three sets are 0.79 ± 0.08, 0.66 ± 0.06, and 0.42 ± 0.06 for 
the variable, AGN, and Belt sources, respectively. The spectral index 
distributions of the 3 sets are quite consistent according to the K and T 
tests, so the slightly softer average index of 2.43 ± 0.05 for the variable 
sources, compared with 2.23 ± 0.02 for the Belt and AGN sources, 
is not significant. Taking into account the detection biases (Grenier, 
2000), the spatial distribution of the variable sources implies an origin 
in a thick Galactic disk with a scale height of 1.3 =L 0.6 kpc. A fit to the 
whole, persistent and variable, source sample with various combinations 
of Galactic distributions nearly equally shares the sources between the 
Gould Belt and a large scale height component, with a 3.4<r and a 4.2a 
improvement over the single distribution fits. At typical distances of 

5 to 10 kpc, the luminosities of the variable sources range from 2 to 
30 x 10 34 erg s 1 sr 1 above 100 MeV. They exhibit large luminosity 
Lry/Lx ratios of a few hundred. 

The compact objects likely to power sources high above the plane 
include ms pulsars and microquasars with a low-mass star companion, 
both having migrated away from the Galactic plane or escaped from 
globular clusters. None of the 7-ray sources coincides with a globular 
cluster. Pulsed 7-rays have been detected from the ms pulsar PSR 
J0218+4232, in phase with the radio and X-ray peaks (Kuiper et al., 
2002). This object shares many traits with the halo sources: a distance 
of 5.7 kpc and an altitude of 1.6 kpc, a luminosity of 1.6 x 10 34 erg s~ x 
sr" 1 and a spectral index of 2.6 above 100 MeV. PSR J0218+4232 does 
not belong to a globular cluster either. Yet, no long-term variability is 
expected from theory (Zhang h Cheng, 2003). In the next section, we 
explore whether low-mass microquasars can produce the halo sources 
despite their intrinsic faintness and softness compared with the young 
high-mass systems. Their stellar and disc intensities are reduced by 4- 

6 and 1 order of magnitude, respectively, and their thermal emissions 
peak a decade or two lower in energy. 



3. Emission from low-mass microquasars 

Microquasars are found up to very high Galactic latitude and height 
above the plane (Mirabel et al., 2001). For instance, XTE J1118+480, 
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at b = 62° and a distance of 1.8 ± 0.6 kpc, lies at 1.6 ± 0.5 kpc above 
the plane (McClintock et al., 2001a). The central object accretes matter 
through a disc from a low-mass star via Roche lobe overflow. McClin- 
tock et al. (2001a) and Wagner et al. (2001) constrained the companion 
spectral type to be between K5 V and Ml V. A large mass function, 
f(M) ss 6 MQ, strongly suggests that the compact object is a black 
hole in a fairly compact binary system with a short orbital of 4.1 hr 
(Wagner et al., 2001). The disc may be precessing under the stellar tidal 
influence (Torres et al., 2002). The rapid and correlated UV-optical-X- 
ray variability in the low-hard state is interpreted as a signature of the 
strong coupling between a hot corona and a jet emitting synchrotron 
radiation up to, at least, the UV band (Hynes et al., 2003; Chaty et 
al., 2003; Malzac et al., 2004). So, pairs in the jet can reach 10 GeV 
in a modest magnetic field of 10 G. Whether the synchrotron emission 
extends to hard X rays is quite possible since the data requires the jet to 
dominate the energy budget and drive the variability, so jet pairs could 
reach hundreds of GeV. The outflow has remained steady through the 
outburst evolution (Chaty et al., 2003). The coronal emission extends 
to ~ 150 keV (McClintock et al., 2001b). The optical to hard-X-ray 
data can me modelled by the Comptonisation in a hot corona or in the 
inner accretion flow of the soft photons emitted by the outer cold disc 
with an inner radius of ~ 55-Rg c h w and a temperature of ~ 24 eV (Esin 
et al, 2001; McClintock et al., 2001b; Malzac et al., 2004). 

XTE J1118+480 therefore serves as a good example for a low-mass 
microquasar and we adopt its characteristics as input to our model (see 
table I). The spectral energy distributions (SED) of the thermal stellar, 
disc and coronal components, taken from the afore mentioned publi- 
cations, are displayed in Figure 3. a. The case of an F star companion 
is also considered in Figure 3.b. We calculate the EC emission from 
the three radiation fields assuming a population of e + — e~ pairs in a 
persistent, cylindrical jet, with a power-law E~ p distribution in number 
density per unit energy between E em \ n and E emax . The jet is assumed 
to be parallel to the disc axis, at an angle $ to the line of sight. It moves 
with a bulk Lorentz factor Tjet and carries a total power Pj e t = qjetMc 2 . 
The coronal emission is assumed to fill a sphere inscribed in the in- 
ner disc radius. The particles are isotropically distributed in the jet 
frame and see the transformed energy densities from the external seed 
photons. The cross-section for both the Thompson and Klein-Nishina 
regimes is used. The Compton losses in the different regions modify 
the injected electron spectrum, introducing a break in the power-law 
from an index p to p + 1 at the energy at which the cooling time equals 
the escape time. This is effective in the disc radiation field. Another 
important ingredient is the absorption from two-photon pair creation 
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Table I. Parameter set used in the model 



black hole mass 




jVibh =0.0 JVLCy 


mass accretion rate 




a'/T Q v If! -8 trr — 1 

ivi — oa _i_u ivit^j yi 


K-M star bolometric luminosity 




Lkm — 4 x 1U erg/s 


F star bolometric luminosity 




L/Y — -L-0 X 1U erg/S 


star temperature 




ftiKM — i cv ana /ci f — 1.0 ev 


star orbital radius 




.L/* = 1. / x 1U cm 


jet/accretion power ratio 




« D / A //" ^.2 in - 3 i n in- 2 
Qjct = ijet/iWC = 1U 10 1U 


corona luminosity 




T 7 o w i n34 / 

ivcor — f.O X 1U CTg/S 


corona outer radius 




JXcor — 1U CHI 


corona photon index (dNx / dE oc 


h x ) 


_, 1 Q 

Clcor = l-O 


corona cut-off energy 




£: cor = 150 keV 


disc luminosity 




^disc = 8.6 x 10 35 erg/s 


disc temperature 




kT disc = 24 eV 


initial jet radius 




Rjet = 1.9 x 10 7 cm 


jet bulk Lorentz factor 




r jct = 3 to 10 


jet viewing angle 




$ = 1° to 30° 


jet electron index (dN e /dE oc E~ 


P ) 


p = 2 to 3 


maximum electron energy 




£ e max = 5 GeV to 5 TeV 


minimum electron energy 




E emin = 1 to 5 MeV 



in the ambient radiation (Romero et al., 2002). It turns out to be quite 
effective near the disc while the coronal and stellar fields are optically 
thin to the 7 rays. Back to the observer's frame, and to observed ener- 
gies shifted by the Doppler factor D = r^(l — /5j e t cos < 1 ) )~ 1 , the stellar 
and coronal contributions to the SED eL e are amplified by D 2+p and 
the disc one by D 2+P (l — cos c£)( 1+ p)/ 2 because of the anisotropic seed 
field (Dermer et al., 1992; Dermer & Schlickeiser, 2002). The former 
amplification factor peaks along the jet axis and the latter near <I> = 15° 
for Tjct ~ 3 and p = 2 to 3. The jet synchrotron emission is only 
amplified by L>( 3 +p)/ 2 . 

This scenario differs from that of Georganopoulos et al. (2002) who 
have imposed a much lower energy cutoff to the electrons so that the 
hard X-ray emission result from their EC interactions with the stellar 
and disc photons rather than from a hot inner accretion flow or from 
a hot coronal plasma energized by magnetic flares above the disc. 

Figure 3 shows the SED obtained per steradian in the laboratory 
frame, adding the contributions from the three external photon fields. 
The coronal component (in the Klein-Nishina regime) is preponderant 
in the COMPTEL 3-30 MeV band whereas the disc component (in the 
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Thomson regime) takes over above 100 MeV. This is why the maximum 
7-ray luminosity is reached for $ close to 15°, reflecting the angle 
dependance of the disc amplification factor. The disc component is 
rather weak because of the local absorption. It has a photon spectral 
index of 2.5 between 0.1 and 1 GeV. The stellar component is small at 
all viewing angles, so the star spectral type has little impact on the 7- 
ray flux. These results show that, even though the spectral index in the 
EGRET band matches that of the unidentified sources, the maximum 
predicted luminosity, L max ~ 4 x 10 29 (i?/100 MeV) -0 - 5 erg s -1 sr" 1 , is 
5 orders of magnitude too faint to account for the halo source flux at 
distances of 5 to 10 kpc. 

The relative contributions of the three EC components change in 
the case of an extreme microblazar where the bulk Lorentz factor Tj e t 
reaches 10, electron energies extend to 5 TeV, and where the jet axis 
is close to the line of sight (<£ = 1°). The generation of such a highly 
relativistic outflow, with rj ct > 10, has been recently observed from 
Circinus X-l, a neutron star with a stellar-mass companion (Fender et 
al., 2004). Particle energies as high as 10 TeV have been inferred in the 
jets of the stellar-mass microquasar XTE J1550-564, 0.1 pc away from 
the black hole, even though the jet is seen deccelerating in the inter- 
stellar medium (Corbel et al., 2002). High-mass extreme microblazars 
are unlikely counterparts to the low-latitude EGRET sources because 
of their extreme brightness. Bosch-Ramon et al. (2004) predict a lumi- 
nosity of 3 x 10 37 erg s" 1 sr _1 at 1 GeV, well in excess of any EGRET 
source. Such an event would appear as bright as Vela at 5 kpc. Figure 
4 shows the result for a low-mass system. The soft EC emission from 
the corona predominates up to several hundred MeV for an F star and 
to several GeV for a K-M star, beyond which the harder, E~ 2 ^, stellar 
component takes over. The disc EC emission is negligible because of the 
lesser amplification at small aspect angle. The 35 times larger energy 
density that the jet encounters around an F star compared with the 
K-M one results in a modest luminosity increase by a factor < 5 in the 
EGRET band. Anyhow, the very low luminosities, as well as the very 
soft spectra, with photon indices ranging from 3 to 4, however, cannot 
account for the EGRET source characteristics. Nor would these systems 
be detected above 100 GeV by the new-generation Cerenkov telescopes 
since the predicted luminosity falls orders of magnitude below their 
sensitivity at a few kpc distance. 

These objects can remain quiescent for years or decades before bright- 
ening by as much as 10 in X rays in a week. One should, however, 
keep in mind that the adopted keV luminosity of ~ 10 35 erg s _1 sr _1 
corresponds to such a flaring state, i.e. to a bright X-ray source with a 
vF v flux of 2.6 x 10~ 4 (D/5kpc)- 2 MeV cm" 2 s _1 . The GeV emission 
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Figure 3. Spectral energy distribution of the EC emission from the jet of a micro- 
quasar with a K-M star (a) or an F star (b) companion, seen at angles of 5° (A), 15° 
(B), and 30° (C) from its axis, for rj ct = 3, <j e t — 1%> an d an E~ 2 3 pair spectrum 
between 1 MeV and 5 GeV. 



being mainly due to the disc or corona (depending on the more or 
less extreme conditions in the jet), the predicted luminosity should 
be considered as an upper limit for comparison with the unidentified 
sources. Therefore, EC emission in low-mass microquasars largely fails 
to explain the variable EGRET sources at large scale heights. 

In contrast with high-mass systems where the external radiation 
energy density largely surpasses the magnetic one, SSC emission is 
likely to dominate even for a modest field strength of ~ 10 G in the jet 
for which the magnetic energy density compares with that of a K-M 
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Figure 4- Spectral energy distribution of the EC emission from the jet of an extreme 
microblazar with an F or K-M star companion, seen at 1° from its axis, for rj ot = 10, 
(fret = 1%, and an E~ 2 3 pair spectrum between 5 MeV and 5 TeV. 
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Figure 5. Spectral energy distribution of the synchrotron and SSC emission from 
the jet of an extreme microblazar, with a 10 G magnetic field, seen at 1° from its 
axis, for Tj et = 10, q^t = 0.1%, and an E~ 2 ' 3 pair spectrum between 5 MeV and 50 
GeV. 



star. Applying the model developed by Bosch-Ramon et al. (2004) to a 
cylindrical jet, one gets luminosities in excess of the necessary EGRET 
ones (see Figure 5). The synchrotron part of the spectrum is 100-1000 
times fainter than the near-IR to UV synchrotron emission recorded 
from XTE J1118+480 during the outburst, the adopted electron index 
in Figure 5 being softer than that of XTE J1118+480 to avoid too hard 
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spectra at high energy. The latter softens beyond 10 GeV because of 
Klein-Nishina effects. One should keep in mind, however, that 77 ab- 
sorption against the disc photons is not included although it efficiently 
limits the emerging 7-ray flux. Adiabatic losses in an expanding jet 
would also modify the result although not severely. Nevertheless, SSC 
emission in a low-mass microquasar appears as an interesting possibility 
to be investigated in greater detail, bearing in mind that modelling the 
high-energy radiation from these fascinating objects is severely limited 
by the highly uncertain choice of the jet bulk motion and magnetic 
field. 
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